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We use both Eilenberger-Usadel and Ginzburg-Landau (GL) theory to calculate the superfluid's 

temperature-dependent kinetic inductance for all currents up to the depairing current in thin and 
narrow superconducting films. The calculations apply to BCS weak-coupling superconductors with 
isotropic gaps and transport mean-free paths much less than the BCS coherence length. The kinetic 
inductance is calculated for the response to a small alternating current when the film is carrying a 
do bias current. In the slow-experiment/fast-relaxation limit, in which the superconducting order 
parameter quasistatically follows the time-dependent current, the kinetic inductance diverges as the 
bias current approaches the depairing value. However, in the fast-experiment/slow-relaxiation limit, 
in which the the superconducting order parameter remains fixed at a value corresponding to the dc 
bias current, the kinetic inductance rises to a finite value at the depairing current. We then use 
time-dependent GL theory to calculate the kinetic impedance of the superfluid, which includes not 
only the kinetic reactance but also the kinetic resistance of the superfluid arising from dissipation 
due to order-parameter relaxation. The kinetic resistance is largest for angular frequencies w obeying 
uJTs > 1, where Ts is the order-parameter relaxation time, and for bias currents close to the depairing 
current. We also include the normal fluid's contribution in deriving an expression for the total kinetic 
impedance. The Appendices contains many details about the temperature-dependent behavior of 
superconductors carrying current up to the depairing value. 

PACS numbers: 74.78.-w,74.78.Na,74.25.F- 



I. INTRODUCTION 

The kinetic inductance, arising chiefly from the ki- 
netic energy of the superfluid, plays an important role in 
superconducting devices fabricated using thin and nar- 
row superconducting films. In such cases the kinetic 
inductance is generally much larger than the geomet- 
ric inductance arising from stored magnetic energy. 
For example, the kinetic inductance plays a prominent 
role in determining the reset time of superconducting 
single-photon detectors (SSPDs) fabricated with mean- 
dering superconducting lines. ''"^ Various calculations of 
the kinetic inductance, relevant to the performance of mi- 
crostrip rcsonators"'^'' and microwave kinetic indiictance 
detectors (MKIDs),^^ have been carried out using (a) the 
London equations neglecting the current-induced sup- 
pression of the order parameter,^"^'^^ (b) the Ginzburg- 
Landau (GL) equations, (c) a GL-inspired London- 
equation approach accounting for the current-induced 
suppression of the order parameter, ^^'^^ and (d) the BCS 
theory.^'^*^ Our goal in this paper is to present theoretical 
calculations of the kinetic inductance for all temperatures 
and for all currents up to the depairing current for sample 
dimensions and properties applicable to present experi- 
mental studies of SSPDs'^^'^^ and micro-resonators.^ Be- 
cause these studies have used thin high-resistance films 
of NbN,i'2>7-9,i8 Nb,2 NbTiN,i9 and TaN2"^2i -^^ ^^i^ ^Yiviy 

limit, we adopt an isotropic s-wave BCS description, al- 
though many of our results can be extended to apply 
under more general assumptions. 

We consider thin {d <C Aq) superconducting films of 
width W much less than the two-dimensional screen- 
ing length (Pearl Icngth^^) A = 2Ao/rf, where Aq is 
the temperature-dependent weak-field London penetra- 



tion depth and d is the film thickness. The condition 
W <^ K guarantees that the self-field generated by the 
current has a negligible effect upon the current density j, 
which therefore flows with the same spatial distribution 
as in the normal state. Moreover, this condition also 
guarantees that the inductance L = L„ -|- is dom- 
inated by the kinetic inductance of the superfluid ife, 
which is typically larger than the geometric inductance 
Lm (associated with the energy stored in the magnetic 
field) by a factor of order K/W.^ We focus on the cal- 
culation of the superfluid's kinetic inductivity Ck- For 
a long strip of length £, width W , and thickness d, the 
kinetic inductance is Lk = Ck^/Wd. 

When the superconductor carries such a low current 
that the superconducting order parameter is not signif- 
icantly suppressed, the electromagnetic behavior is well 
described by the London equation, and the kinetic energy 
density of the superfluid can be expressed as^ 

rr 1 2 1 / \ -2 1 \2 -2 1 /- -2 

(1) 

where rigo is the superfluid density, m the electron mass, 
Vs the superfluid velocity component in the x direction, 
js = —UsQeVs the supercurrent density component in the 
X direction, and — e the electron charge, such that the 
kinetic inductivity of the superfluid obeys^^'^^ 

Uq (T) = Mo Ag (T) = m/n^oe^ • (2) 

The subscripts on rigo, Aq, and £fco are a reminder 
that these quantities apply in the limit as js — ^ 0. The 
simple relationship given in Eq. (2) has made possible 
the determination of Ao(r) vs T in YBa2Cu307_5 from 
kinetic- inductance measurements.^® 
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When the superconductor carries high currents, how- 
ever, calculation of the superfluid's kinetic inductance 
becomes more complicated, especially when the current 
density approaches the depairing value ja, which we cal- 
culate in Sec. II using the mean-field Eilenberger-Usadel 
theory.^^'^^ In high currents it is no longer possible to de- 
fine the kinetic inductance by considering only the stored 
kinetic energy density as in Eq. (1), because, as exam- 
ined in detail in Appendix D, increasing js to large values 
also affects the superconducting condensation energy by 
suppressing the superconducting order parameter. 

An additional complication is that n,, which depends 
on the superconducting order parameter, can change only 
on a time scale slower than a characteristic (but compli- 
cated to determine^'*'^^) relaxation time Tg. As a conse- 
quence, Hs may or may not be able to follow the changes 
in js and Vg . 

Calculations of the current dependence of jC-k sue sim- 
plified in two limiting cases: (a) slow experiments (fast 
relaxation, Sec. Ill), in which jg and Vg vary on an exper- 
imental time scale Texp much longer than the relaxation 
time Tg, such that the order parameter and the superfluid 
density quasistatically follow js and Vs, and (b) fast 
experiments (slow relaxation, Sec. IV), in which js and 
Vs change so rapidly about their time averages js and Vs 
(on a time scale r^^p much shorter than t^) that the or- 
der parameter and Us cannot track the time dependence, 
and Hs remains very close to the value corresponding to 
js and Vs- To provide an approximation to the transi- 
tion between these two limiting cases, in Sec. V we em- 
ploy a simplified phenomenological model based on the 
time-dependent GL (TDGL) equations^'' to calculate the 
complex impedivity due to the superfluid. In Sec. VI we 
include the normal-fluid's resistive contribution to the to- 
tal complex impedivity, and in Sec. VII we provide a brief 
summary and discussion of our results. Various details 
of the calculation are included in Appendices A-E. 



II. SUPERFLUID- VELOCITY DEPENDENCE 
OF THE SUPERCURRENT DENSITY AND THE 
DEPAIRING CURRENT DENSITY 

We next study the behavior of the order parameter 
and the current- carrying capacity of superconductor in 
response to a slowly varying applied current. Since we 
are interested in a quasistatic theory valid in supercon- 
ductors with a short normal-state mean-free path and at 
all temperatures in the superconducting state, we employ 
the quasiclassical Eilenberger^'^ theory as formulated by 
UsadcP^ for the dirty limit. Consider a superconducting 
strip extending along the x direction when the current 
is uniform. Let js and As = mVs/e denote the x com- 
ponents of the supcrcurrent density js and the gauge- 
invariant vector potential Ag = A + (0o/27r)V7, where 
A is the gauge-dependent vector potential, (f)o = h/2e 
the superconducting flux quantum, and 7 the gauge- 
dependent phase of the superconducting order param- 



eter. 

The supercurrent density can always be expressed as 
js = —UseVs, but in general is a function of the super- 
fluid velocity Vs and has the value Usq when Vs = but 
decreases monotonically to zero as \vs\ increases. For 
positive js and negative values of Vs, the supercurrent 
density js = nse\vs\ initially increases linearly as a func- 
tion of \vs\, reaches a maximum jd{T) (the depairing or 
pair-breaking current density) at \vs\ = Vd{T), then de- 
creases to zero at \vs\ = Vm{T). 

When the superconductor is current-biased, only the 
portion of the curve js vs \vs\ for < \vs\ < Vd{T) is 
accessible. On the other hand, following a suggestion 
by Fulde and Ferrell,^° Bhatnagar and Stern^^'^^ showed 
that it is possible to probe experimentally the shape of js 
vs \vs\ even for Vd{T) < \vs\ < Vm{T) using a multiply- 
connected sample geometry. In this paper we first exam- 
ine the behavior of js over the full range of values of Vg, 
but later in applying these results to study the kinetic 
inductance we limit our attention to the current-biased 
case in which js is a single- valued function of Vs in the 
range Q < js < jd- 



A. Depairing current density calculated from the 
Usadel equations 

For the problem at hand the Usadel equations can be 
written as^^ 

-hD{GF' - FG')' = 2AG - 2hWnF , (3) 
G^ + \Ff = l, (4) 

Ai„^ = 2,..r|(^-F), (5) 

00 

js = -ATreN{0)DkBT ^ lmF*F' , (6) 

71=0 

where A is the superconducting order parameter, D ~ 
VpT/3 = Ui?£/3 is the diffusivity, vp is the average ve- 
locity of electrons at the Fermi surface, r is the normal- 
state transport lifetime, £ is the mean-free path, hujn = 
(2n -|- VjiikBT is the Matsubara frequency, N[Q) is the 
density of Bloch states of one spin at the Fermi level, T 
is the temperature, and Tco is the zero-current transition 
temperature. The primes in Eq. (3) denote differentia- 
tion with respect to x. These equations describe super- 
current flow in a superconductor with an s-wave isotropic 
gap in the weak-coupling limit of the BCS theory. 

However, this mean-field theory does not account for 
the possibility that one- or two-dimensional fluctuations 
could grow to produce phase slips or vortex crossings. 

Since W ^ A, we can neglect the self-field of the cur- 
rent and choose a gauge for which we may replace the 
gauge- invariant gradient V -|- 2-jTiA/(po by xd^- Looking 
for solutions of the form A = Age*"^^, F = F„qe*'^, and 
G = Gnq, where q is the gradient of the phase of the 
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order parameter, we find that Eqs. (3) and (4) become 



G 



nq 



f"^ = 1 

^ nq ^ ■ 



(7) 
(8) 



where Q = TiDq^ /2. 

As noted by Maki,^^ these equations are equivalent to 
those of the Abrikosov-Gor'kov (AG) theory^^ for pair- 
breaking scattering, except for the replacement of the 
AG spin- flip scattering rate l/rm by Dq^ /2. Our re- 
sults therefore share many properties with those of the 
AG theory. For example, we show that the transition 
temperature depends upon q and decreases mono- 
tonically from its value Tco at g = to zero at a critical 
value of q given by g™(0) = 1/^(0) = {^iot/^)-^'^ , where 
^0 = Tivp /tt^o{0) is the BGS coherence length.^^ For a 
fixed value of q, the order parameter Ag(T) is nonzero 
only for temperatures T less than T^q] equivalently, for 
a fixed temperature T, the order parameter Aq(T) is 
nonzero only for values of q less than qm{T) = 1/^(T). 
Note that since the phase of the order parameter is 
7 = qx, the gauge-invariant vector potential is related to 
q and the superfluid velocity via Ag = {(f)Q/2iT)q = mVs/e. 

Introducing Unq = Gnq/Fnq, we find that Eqs. (7) and 
(8) can be written as 
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FIG. 1: [A,(r)/Ao(0)]2 vs q/qm{0) obtained from Eqs. (11) 
and (A2) for various values of t = T/Tcq. The dotted lines 
show the linear behavior as [Aq(T)/Ao(0)]^ — ^ in the limit 
as g qm{T). 
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and 



^ nq — 



2 

'nq 



(9) 



(10) 



where rj = n + 1/2, e = Ag/27rfcBT and C = Q/^q- Unq 
(which depends implicity upon T) can be obtained for 
arbitrary values of rj, e, and ( by solving Eq. (9) as a 
quartic equation [see Appendix A]. 

With the introduction of Unq, two equations remain 
to be solved. The self-consistency equation (5) in the 
presence of the current becomes, for general values of 
Ag, CO, Q, and T, 



In- 



t 



^ I n+1/2 

71=0 » ' 



2 

nq J 



In general, Ag(T) must be obtained by numerically solv- 
ing Eq. (11) using Eq. (A2), but the results can be 
checked against analytic results obtainable in the limits 
of g ^ and g ^ qra{T). Figure 1 shows [A,(r)/Ao(0)]2 
as a function of q/qm{0) for a series of values of the re- 
duced temperature t = T/Tco- The sums for t > 0.1 were 
evaluated by summing n from to 500, but for t = we 
used the analytic results in Eqs. (23)-(25). 

From the current equation (6) we find that when the 
superfluid velocity Vg is in the x direction, the supercur- 
rent density in that direction is^^ 



3sq{T) = -nsq{T)eVs. 



(12) 



From Eqs. (6), (10) and (12) we obtain a general expres- 
sion for the g-dependent superfluid density. 



_ 87TmN{0)DkBT ^ 1 

n=0 



(13) 



When 5 ^ 0, we have u„o = 27TkBT{n + l/2)/ Ao{T) [see 
Eq. (9)], and when this is used in Eq. (13), evaluation of 
the sum yields 



nso{T) = 
such that 



27rm7V(0)£)Ao(T) 



h 



«so(0) = 



tanh 



V 2kRT J'' 



2kBT . 
27rmiV(0)DAo(0) 

n 



(11) and 



nsq{T) ^ AkuT 
n,o(0) Ao(0) 



o(0) ^ 1 + ti 



(14) 



(15) 



(16) 



In general, nsq{T)/nso{Q) must be obtained by numer- 
ically solving Eqs. (11) and (16) using Eq. (A2). Figure 
2 shows nsq{T)/nso{Q) as a function of q/qm{0) for a se- 
ries of values of the reduced temperature t = T/Tcq. The 
sums for t > 0.1 were evaluated by summing n from 
to 500, but for t = we used the analytic results in Eqs. 
(26)-(27). 

As shown in Fig. 2, nsq{T)/nso{0) depends upon q and 

vanishes at g = qm{T) (see Appendix C). The corre- 
sponding g-dependent penetration depth \{T) can be 
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FIG. 2: ■n„,(T)/n,„(0) = \l{Q)/\l{T) vs g/g™(0) obtained 
from Eqs. (11), (16), and (A2) for various values olt = T/Tco. 
nsg{T)/nso{Q) as q ^ qm{T). 



obtained from 



^.o(0) A2(T)- 



(17) 



Note, however, that current-biased experiments can ac- 
cess values of q only up to qd{T), where the magnitude 

of the current density reaches the depairing Hmit jd(T). 
The general expression for the supercurrent density is 



A.. 



m 



(18) 



Prom Eq. (12) or (18) we see that, because jsq{T) is 
the product of nsq{T) (a monotonically decreasing hmc- 
tion of q) and eVs = e^As/m = ehq/2m, the magnitude 
of jsq{T) reaches a maximum, called the depairing cur- 
rent density jd{T), when q = qd{T), where < qd(T) < 
qm{T). We define jq{T) as the magnitude of jsq{T) nor- 
malized to nso{Q)eVm{^) = ?!>o/27r/ioAo(0)^^(0), such that 



n, 



AT) 



risoiO) gm(0)' 



(19) 



The maximum value of jq(T) vs q is the normalized de- 
pairing current density jd{T). 

In general, jq{T) must be obtained numerically from 
Eqs. (11), (16), (19), and (A2), but the results can be 
checked against analytic results obtainable in the limits 
f — >■ and f — >■ 1, to be discussed in more detail later in 
Sees. 11 B and 11 C. Figure 3 shows the general behavior 
of jq{T) as a function of q/qm{^) for a series of values 
of the reduced temperature t = T/Tco- The points la- 
bel the values of jg and q corresponding to the depair- 
ing current density jd and qd- The solid curve in Fig. 4 
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FIG. 3: jq{T) vs q/qm{0) obtained from Eqs. (11), (16), (19), 
and (A2) for various valus oi t = T/Tco- The points label 
the values of jq and q corresponding to the depairing current 
density jd and qd. Current-biased experiments probe only the 
portions of the curves to the left of these points. jq{T) — ^ 
as q^ qm{T). 
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FIG. 4: b'd(T')/id(0)]^/^ (solid) vs t = T/Tco obtained 
numerically from Eqs. (11), (16), (19), and (A2). The 
dotted line shows the behavior of jd{T) in the GL limit 
near Tco [Eq. (41)], and the long-dashed curve shows the 



approximation' 



j4t)/M0) ^ (1 - n 



2\3/2 



The short- 



dashed curve shows the variation of Pd{T) [Eq. (20)] from 
0.475 at t = to 0.385 at t = 1, and the dot-dashed curve 
shows the variation of kd{T) [Eq. (21)] from 0.595 at t = to 
0.544 at t = 1. 



shows [jd(T)/id(0)]2/3 = [jrf(r)/3d(0)]2/3 as a function 
oi t = T/Tco, and the dotted line illustrates how jc{T) 
approaches the (1 — t)^/^ behavior in the GL regime close 
to Tco- 

For all temperatures below T^q, estimates oi jd{T) can 
be obtained from 



Jd{T)=Pd{T) 



27rMoAg(T)^(r)' 



(20) 



where Pd{T) is a dimensionless function defined by Eq. 
(20). The dashed curve in Fig. 4 shows Pd{T), which 
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is obtained numerically from all the other quantities in 
Eq. (20), varies from 0.475 at T = [Eq. (30)] to 0.385 
as T Teo [Eq. (42)] with a maximum of 0.483 at t = 
T/Tm = 0.17. 

Similarly, estimates of the depairing current density 
also can be obtained for all temperatures from 



jd{T) = kd{T)H,{T)/Xo{T), 



(21) 



where the dimensionless quantity kd{T), defined by Eq. 
(21) and shown by the dot-dashed curve in Fig. 4, varies 
from 0.595 at T = [Eq. (30)] to 0.544 as T -i^ [Eq. 
(42)] with a maximum of 0.608 &tt = T jT^a = 0.21. The 
values of fc^ shown in Fig. 4 were obtained from 

t.m=0.5.5¥j2fffl^ (22) 
3<0) ^c{T) Ao(0) 

via Eqs. (32), (E2) and (B2), where Hc{T) is the 
temperature-dependent bulk thermodynamic critical 
field (see Appendix E). 



B. Depairing current density at zero temperature 

At T = 0, the q dependence of Ag(0) can be obtained 
by converting the sum in Eq. (11) to an integral over u„g. 
The result is^^'^^ 



A,(0) 
Ao(0) 



exp(-7rCo/4), < Co < 1, (23) 

exp[-(Co sin-i Co"' - \/l-Co"')/2 
-cosh-^Co], Co>l, (24) 



where 



Co 



hPq^ ^ 1/ q \2 Ao(0) 
2Ag(0) 2\qm{0)J Ag(0)- 



(25) 



Figure 5 shows [Aq(0)/Ao(0)]^, obtained from numerical 
solution of Eqs. (23)-(25), as a function of q/qm{0) = 

The q dependence of nsg(O) can be obtained in a similar 
way. The result is^^ 



= exp(-7rCo/4)(l - 4Co/37r), < Co < 1, (26) 



= exp(-^Co/4){:— :2 [(1 + 2Co)\/Co' " 1 " 2Co 



SttCo^ 

+ -sin-iCo"'}, Co>l, 

TT 



(27) 



where Co is given by Eq. (25). nsq{0) is shown as the 

dashed curve in Fig. 6. The g-dcpendcnt penetration 
depth at zero temperature Xq{0) can be obtained from 



n.o(0) A2(0)" 



(28) 
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FIG. 5: [Aq(0)/Ao(0)]^ vs g/gm(0) obtained from Eqs. (23)- 
(25). The dashed line shows the linear behavior 6[1 — ^/^^(O)] 
as g ^ 9m(0). 
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FIG. 6: Reduced superfluid density nsq(0)/nso(0) = 
Ao(0)/Aq(0) (dashed curve) vs q/qm{0) — Vs/vm{0), obtained 
from Eqs. (26) and (27), and normalized q-dcpcndcnt super- 
current density jq{0) (solid curve) vs q/qm{0) — Vs/vm{0), 
obtained from Eq. (29). The black point and dotted lines 
indicate the maximum jq{0), the normalized depairing super- 
current density jd = 0.475, which occurs at qd/qm{0) = 0.689. 



To obtain the depairing current density, consider the 
g-dependent (or Wg-dcpendent) supercurrent density at 
T = [Eq. (12)], normalized to -nso(0)ew„(0) = 



(0) 



«so(0) Vg„(0) 



(o^foO' 



(29) 



shown as the solid ciirve in Fig. 6. The point and dot- 
ted lines show the maximum .7^(0), the normalized de- 
pairing supercurrent density jd(0) = 0.475, which occurs 
at c/rf(0)/g™(0) = 0.689 and Co = 0.300. The resulting 
zero-temperature depairing supercurrent density can be 
expressed in several ways: 



Jd(0) 



0.475- 



'27r/.oA2(0)C(0)' 

1.491iV(0)e[Ao(0)]3/2/D7?^, 
0.595i?c(0)/Ao(0). 



(30) 

(31) 
(32) 
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Equation (31) coincides with the rcsuh given for the de- 
pairing supercurrent density in Ref. 39. Here Hc{T) is 
the temperature-dependent bulk thermodynamic critical 
field (see Appendix E), and iJc(O) = Ao(0)^A^(0)//io- 



C. Depairing current density in the GL regime 

As shown by Gor'kov^*^, the GL theory of 
superconductivity^^ is derivable from the micro- 
scopic BCS theory^^ at temperatures T very close to 
the transition temperature Tco- This assures us that we 
also can apply the Usadel equations to recover the well- 
known GL results. For T close to T^q or, equivalently, 
for q close to qrn{T), the g-dependent order parameter 
Aq (T) becomes small, and it is useful to expand Unq and 

+ ulg in powers of e = Aq{T) /2TrkBT: 



(r? + a)3 {tj + a)4 



where 



a = 



T] + a 2 
+0(e^), 

e-^qljT) ^ O.UOqliT) 
At g^(0) t ql{Q)- 



(33) 



(34) 



Substituting this into Eq. (11) and keeping only 
the lowest order terms, since we know that e = 
Ag(r)/27rfcsT < 1 when 1 - t < 1, where t = T/T^, we 
obtain 



In 



oo ^ 



1 



+ 



(n+l/2 + Q)3 (n+l/2+a)4j 2 



(35) 



The sums can be expressed in terms of digamma func- 
tions and their derivatives. When l — t"^ 1, Eq. (35) has 
solutions only for a ^ 1 and can be expanded as 



l-t = 



2 - 



7C(3)e2 



89^(0) 



(36) 



Solving for e^, dividing by Eq = 2(1 — i)/7C(3), and mak- 
ing use of g^(r) = [8e'^(l - f)/7r2]g^(0) (see Appendix 
C), we obtain 



A|(T) 
Al{T) 



= 1 - 



(37) 



where qm{T) = l/e(T). 

Substituting the expansion of Eq. (33) into Eq. (16), 
keeping only the lowest order terms, we obtain 



nsqjT) ^ \1{T) ^ AliT) 
n,o(r) A2(T) Al{T) 



(38) 



and 



n,o(0) A2(T) 7C(3) V 



,(0) A2(T) 7C(3)V^ 9?„(T)/^ ^^^^ 
The reduced q-dependent supercurrent density becomes 



47re^ g„(T) 



^"^^^ 7C(3)Z(0)(^ qUT)^qm{T) 



(1-i), (40) 



whose maximum occurs at qd{T)/qm{T) = l/-\/3, such 
that (see Appendix C) the reduced depairing current den- 
sity is 

UT) = ^^^(1 - = 1-230(1 - tf'. (41) 

Thus, in the GL regime the depairing current density can 
be expressed as: 



3d{T) = 0.385; 



(42) 



-27rMoA§(T)e(r)' 
where 0.385 =2/3\/3, 

jrf(T) = 3.865Af(0)e[Ao(0)]3/2/D7n(l - tf^^ (43) 

or, since \f2Hc = (j)Q/2-KHo\o^ in the GL theory, 

irf(T) = 0.544ife(T)Ao(T), (44) 

where 0.544 = (2/3)^/2. 

To simpUfy calculations in the GL limit later in Sees. 

Ill and IV we introduce the parameter (f) {Q < (j) < 7r/2), 
such that 

\3sq\/jd = sm4>, (45) 

q/qm{T) = (2/y3)sin(<^/3), (46) 

f = [l-|-2cos(2(/)/3)]/3. (47) 



III. SUPERFLUID KINETIC INDUCTIVITY IN 
SLOW EXPERIMENTS (FAST RELAXATION) 

We now can calculate the kinetic inductivity of the 
superfluid measured in slow, low-frequency (or, equiv- 
alently, fast relaxation) current-biased experiments, in 
which both j, < u{T)] and v, [\v,\ < Vd{T)] 

vary on a time scale Texp much longer than the GL 
relaxation time Ts-^^ In a one-dimensional conductor 
carrying a uniform current the gauge-invariant electric 
potential P = $ — {(j)Q/2iT)d'y/dt is zero,''^ and the 
electric field along the conductor is = —dAg/dt = 
Ck{q,T)djsq{T)/dt.^^ Since from Eq. (18) we have As = 
(0o/27r)g = —jsq{T)fioX'^{T), taking the time derivative 
and using df/dt = {df /dq)dq/dt we obtain the kinetic 
inductivity of the superfluid for slow experiments. 





-1 


djsq{T) 


[dq\Xl{T)J\ 




dAs 



= «A5(r,f.(iM.) 



(48) 
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where the slow-expcrimcnt function Fs is simply 
Ck{q, T) / ^qXq{T) but expressed as a function of the nor- 
mahzed current density \js\/jd{T) rather than as a func- 
tion of q. While g is a convenient theoretical variable, js 
is a more convenient variable for the current-biased case. 
For l^l < qd{T) and \js\ < jd{T), jg is a single- valued 
function of q, shown in Figs. 3 and 6. 

In the limit of small currents, when q ^ 0, Ck{q,T) 
reduces to Ck{0,T) ~ Cko{T) [Eq. (2)]. However, as can 
be seen from Figs. 3 and 6, \djsq/dq\ decreases mono- 
tonically for increasing values of q and becomes zero 
at the depairing value. Accordingly, as q increases, 
Ck{q,T) starts from Cko{T), increases monotonically, 
and diverges at g = qd{T), where \js\ = jd{T). Because 
Ck{q,T)/Cko{T) diverges as -> jd-. we show in Fig. 7 
the typical dependence of the inverse, CkoiT)/ Ck{q, T) = 
^/Fs{\3s\/3d{T)), ys |is|/jd(T). This figure was obtained 
by (a) evaluating jq{T) [Eq. 19)] and djq{T)/dq numer- 
ically for t = T/Tco = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8, and 0.9, and analytically [see Eq. (40)] in the GL 
limit f — >■ 1, (b) calculating 




0.4 0.6 

X = \.is\/.id 



FIG. 7: Cko{T)/Ck{q,T) = l/F,(|i,|/jd(T)) for slow exper- 
iments vs a; = \js\/jd{T) at t = T/Tco = (solid), t = 0.3 
(dotted) and t ^- 1 (dashed). Note that jCko{T) = Ck{Q,T) = 
MoAg(r). 



Cko{T) n,o(0) 



Ck{q,T) n,o{T) 



9m (0) 



dm 

dq 



(49) and 6 = sm ^ x. For small values of x, 



using Eq. (14) to evaluate nsQ{0)/nso{T), and (c) mak- 
ing a parametric plot of Cka{T)/Ck{q,T) vs \js\/jd = 
jq{T)/jd{T). As shown by the solid curve for t = 0, 
the dotted curve for t = 0.3, and the dashed curve for 
t ^ 1, the behavior of jCko{T)/Ck{q,T) vs \js\/jdiT) 
is not monotonic as the temperature changes, but the 
curves for all other temperatures (not shown) lie in a 
narrow band between the dotted and dashed curves. As 
\.is\/.jd 1, all the curves have an inverse-square-root 
dependence close to that in the GL limit f — >■ 1, 

Cg^{T)/jC^^{q,T) = {2VQ/S)[1 - \js\/jd{T)]'/' ■ (50) 

The curves shown in Fig. 7 can be represented by 
ysn{x) = (1 — x")-^/"' (not shown in Fig. 7), which 
fits the calculated values of y = £ko{T)/Ck{q,T) vs 
X = \js\/jd{T) for < X < 0.97 with 1% accuracy for 
{n,t) = (2.21, 0), (2.21, 0.1), (2.27, 0.2), (2.30, 0.3), (2.28, 
0.4), (2.25, 0.5), (2.22, 0.6), (2.18, 0.7), (2.16, 0.8), (2.13, 
0.9), and (2.11, t 1). 

To calculate the kinetic inductivity in the GL limit 
shown by the dashed curve in Fig. 7, it is convenient 
to use the parametric relations x = \js\/jd = sine/) and 
y = Cg^iT)/C'^^iq,T) = 2 cos(20/3) - 1, where < (j) < 
tt/2 [see Eqs. (45)-(47)]. The slow-experiment kinetic 
inductivity of the superfluid in the GL limit is 



{q,T)=t,,\l{T)F^'^[ 



\3s 



3d{T) 



where 



{x) = 



1 



2cos(20/3) - 1 



(51) 



(52) 



Ff^{x) = l + lx' + ^^x^ + 0{x% 



(53) 



and F^'"{x) diverges at a; = 1, as noted in Ref. 13. (See 
also the upper solid curve in Fig. 10.) 



IV. SUPERFLUID KINETIC INDUCTIVITY IN 
FAST EXPERIMENTS (SLOW RELAXATION) 

We next consider fast experiments (or, equivalently, 
slow relaxation), in which the current density js{t) = js + 
jsi (t) as a function of the time t changes rapidly about its 
time average js on a time scale Tgxp much shorter than 
the relaxation time Tg.^^ In this case neither the order 
parameter Ag nor the g-dcpcndent penetration depth 
can follow the time dependence of the current, but in- 
stead they remain frozen to their values at q = q given 
by js = jsq = -As/^i.q\% where As = mVs/e = nq/2e. 
From js{t) = -As{t)/iJLoXl, As{t) = A^ + ^i(t), and 
E{t) = -dAs{t)/dt = Ck{q,T)djs{t)/dt, we find that the 
kinetic inductivity of the superfluid in fast experiments 
is 



Ck{q,T) = MoA?(T) = ^oA^(T)F/ (^) , (54) 

which can be evaluated mimcrically using Eq. (16) as 



Cko{T) nsq{T) 



n=0 



"0 „=o ^ + ""9 



(55) 



(See Fig. 2.) Here the fast-experiment function Fj is 
simply Lk{q,T) / iiq\q{T) but expressed as a function of 
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the normalized current density \js\/jd{T) rather than 
as a function of q. When q — > 0, Ck{q,T) reduces to 
£fc(0,T)=£feo(r) [Eq. (2)]. 

Shown in Fig. 8 is the typical dependence of 
Ck{q,T)/Cko{T) vsjjs\/jd{T). This figure was obtained 
by (a) evaluating jq{T) and Ck{q,T)/ CkoiT) numeri- 
cally for t = T/Tco = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8, and 0.9, and analytically in the GL limit f — )• 1 
[see Eqs. (45) and (47)], and (b) making a parametric 
plot of Ck{q,T)/Cko{T) vs \j,\/jd = jg{T)/jd{T). As 
shown by the solid curve for f = 0, the dotted curve for 
t = 0.3, and the dashed curve for t 1, the behav- 
ior of £/c(g, T)/£fco(r) vs \js\/jd is not monotonic as the 
temperature changes, but the curves for all other tem- 
peratures (not shown) lie in a narrow band between the 
dotted and dashed curves. As \js\/jd — >■ 1, all the curves 
approach their limiting values in the range 1.41 - 1.50 
(solid symbols in Fig. 8) with infinite slope. 

The curves shown in Fig. 8 can be represented by 
yfn{x) = yo — (yo — 1)(1 — a;")-'-/" (not shown in Fig. 8), 
which fits the calculated values oi y = £k{q,T)/£ko{T) 
ys X = \js\/jd, where yo is the value of Ck{q,T) / Cko{T) 
at X = 1, for < a; < 0.97 within 0.5% for {yo,n,t) 
= (1.451, 2.48, 0), (1.448, 2.47, 0.1), (1.422, 2.45, 0.2), 
(1.412, 2.46, 0.3), (1.417, 2.50, 0.4), (1.432, 2.50, 0.5), 
(1.448, 2.50, 0.6), (1.463, 2.50, 0.7), (1.477, 2.50, 0.8), 
(1.490, 2.50, 0.9), and (1.500, 2.50, t 1). 

To calculate Ck{q,T) in the GL limit shown by the 
dashed curve in Fig. 8, we used Eqs. (45)-(47) The kinetic 
inductivity of the superfluid for fast experiments in the 
GL limit is 



£r(,-T) = MoA^(T)Ff^(^), 

^ 3d ■' 



where 



^'^'^^''^ ~ lHx)~ r+2( 



(56) 



(57) 



! cos(2(?i/3) 

and (j) = sin~^ x. As noted in Ref. 13, for small values of 

X, 

Ff\x) = i + ^x' + + 0{x% (58) 

and F'^^{x) approaches 3/2 with infinite slope as a; — )• 1. 
(See also the lower solid curve in Fig. 10.) 



V. KINETIC IMPEDIVITY OF THE 
SUPERFLUID Zus 

In the above sections we discussed the situations when 
Texp/Ts is largo or small. To describe the transition 
between these two limits, for the moment wo restrict 
our attention to the GL regime and employ a simpli- 
fied phenomenological model assuming that the time de- 
pendence of / is determined by simplest version of the 
time-dependent GL (TDGL) equation, 

'2 f (59) 




0.4 0.6 

X = UslUd 



1.0 



FIG. 8: Ck{q.T)JCMT) = F ; (\j s\ / j d{T)) for fast experi- 
ments vs X = \js\/jd(T) ai t = T/TcQ = (solid curve 
and filled circle), t = 0.3 (dotted curve and filled square) 
and t — >■ 1 (dashed curve and filled triangle). Note that 
jCko{T) = noXliT). 



where A'^ = As/{(j)o/2nS^). In the GL regime, the su- 
percurrent density [Eq. (18)] becomes in dimensionless 
quantities 



(60) 



where j'^ = js/(0o/27r/xo^Ag). We now consider experi- 
ments in which the supercurrent density {t) as a func- 
tion of the time t changes about its time average on a time 
scale Texp comparable with the relaxation time Tg . In par- 
ticular, we consider the linear response of the supercon- 
ducting strip to a time-dependent supercurrent density 
given by j'^{t) = i^o where J'^q, the bias current 

current density, is fixed to be in the range < jj^ol < J^, 
and j'^i, the amplitude of the ac current density, obeys 
j'si ^ fso- III this section we assume that the frequen- 
cies are sufficiently low that normal-fluid currents are not 
excited such that the current is all supercurrent. To an- 
alyze the linear response of the reduced order parame- 
ter to the ac current, we substitute = j'^Q + j^^e 
and / = /o + fie"^' (|/i| < /o) into Eq. (59), where 
/o (l/\/3 < /o < 1) is the solution of Eq. (59) in the 
time-independent case when j'^ = j'^Q . We then linearize 
Eq. (59) by neglecting terms of order j'^^ and ff. The 
solution is 



iijt 



/3(6/2-4 + ia;r«)- 



(61) 



T,df/dt = f-r- A'ff. 



From E = -dAJdt, Eq. (60), and j% = /^(l - Q we 
obtain the electric field in the linear-response approxima- 
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tion: 



E 



2/0^ 



^/o(6/o -4 + «^T,)y dt 

Zksjsl = {T^ks + i^ks)jal 
{nks+il^Ck)3sl- 



(62) 

(63) 
(64) 



The coefficient of djsi/dt on the right-hand side of 
Eq. (62) reduces to the slow-experiment inductivity 
£^^{q,T) [Eq. (51)] in the hmit ujTs and to the 
fast-experiment inductivity £^^(g, T) [Eq. (56)] in the 
limit iOTs — >■ 00. 

The complex kinetic impedivity (specific impedance or 
complex resistivity) of the superfluid Zks , here evaluated 
in the GL regime, can be conveniently expressed in terms 
of F^^ [Eq. (52)] and F^^ [Eq. (57)] as 



where 



reff = i^rr«/2. 



(65) 



(66) 



(For a long strip of length £, width W, and thickness d, 
the complex kinetic impedance is Zks = Zks^/Wd.) 

The real part Tiks of the superfluid kinetic impedivity 
is the frequency-dependent resistivity of the superfluid 
due to order-parameter relaxation. Using the the para- 
metric relation x = \ jso\/jd = sin^ as above, we obtain 



n. 



ks 



MoA| /|i,,o| 

G .UJT, 



G{x,UJTs) 

a(<A) 

m 



Ts ^ 3d 

a2(<^) + (a;T«)2' 
4cos(2(^/3) - 2, 

16sin2(0/3) 
l-h2cos(2</)/3)' 



(67) 

(68) 
(69) 
(70) 



G{x, (jjTs) is shown in Fig. 9 as a function of x for various 
values of wr^. In limiting cases, we have 



G{x,UJTs) 



16 2 
27^^ 


64 


243 


[4 2 
[27" 


16 


^81 



a;4 + 0(x^), uJTs > 1, (71) 
a;4 + 0(a;6)l (a;r«)^ 



i^Ts < 1, (72) 

and G{x,ijOTs) approaches 2 with infinite slope as a; ^ 1. 

Although the superconducting strip has zero dc elec- 
trical resistivity, under ac conditions order-parameter 
relaxation contributes to dissipation of energy in a 
manner similar to the way it contributes to flux-flow 
dissipation. 24'^''"'*'' The time-averaged rate of energy dis- 
sipation per unit volume via order-parameter relaxation 
is {l/2)^lksj1l^ which also can be calculated using the 
dissipation function discussed in Refs. 45 and 46. 
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FIG. 9: G{x,ujTs), which describes the alternating-current 
resistivity of the superfluid due to order-parameter relaxation 
[Eqs. (67) and (68)], vs a; = \jso\/jd for ujts = 0.1 (dotted), 
u!Ts = 0.3 (dot-dashed), ujts = 1 (dashed), lots = 3 (long 
dash), and loTs = 00. G{1,uiTs) = 2. 
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FIG. 10: H{x,ujTe), which describes the superfluid's kinetic 
inductivity [Eqs. (73) and (74)], vs a; = Ijsol/jd for ujTs = 
[upper solid curve, for which H{x,0) = F^^{x), Eq. (52)], 
LOTs = 0.3 (dot-dashed), lots — 1 (dashed), cjt., = 3 (long 
dash), and ljts = 00 [lower solid curve, for which H{x,0) = 
Ff^(x), Eq. (57)]. For nonzero lots, H{1,uits) = 1.5 (black 
point). 



The supcirfluid's kinetic; rc;activity is Xks = ^^kj where 
the superfluid's kinetic inductivity £fc is 



Ck =MoA^F(^,wr,), (73) 

16a((?i) sin^(^/3) 



iJ(a;,wTs) 



3+ 



l-F2cos(2(?i/3) L o?{4>) ^ (yjT^y 



(74) 



Shown in Fig. 10 are plots of H{x,ujTs) vs x for various 
values of ojts- As expected, H{x,uits) approaches the 
slow-experiment result F^^{x) when ojTs and the 
fast-experiment result F^^{x) as uoTs — > 00. 

In the above we have used the relatively simple formal- 
ism of the time-dependent GL theory. A reasonable start- 
ing point for a phcnomenological theory of the complex 
kinetic impedivity of the superfluid Zks at lower temper- 
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atures outside the GL regime is to replace the quantities 
F^^ and Ff'^ in Eqs. (65) and (66) by the more gen- 
eral expressions Fg and Ff given in Eqs. (48) and (54). 
Note from Figs. 7 and 8 that Fs = Ck{q,T) / iif)\l{T) 
and Ff = Ck{q,T) / ijlq\q{T) as functions of \is\/id and 
\js\/jd do not differ greatly from their GL counterparts, 
F^^ and F^^. However, as discussed in Ref. 29, there 
are important open questions of what is the meaning of 
the relaxation time Tg and how does one calculate its 
numerical value. 



VI. KINETIC IMPEDANCE INCLUDING THE 
NORMAL-FLUID RESPONSE 

The time-dependent GL theory^^'^^"'*^ is a two-fluid 
theory [generally applicable for frequencies uj well below 
the superconducting gap frequency 2Ag(T)]. In our situ- 
ation the ac current density ji is the sum of the supercur- 
rent density jsi and the normal-fluid current dcmsity j„i, 
the flow of thermally excited quasiparticles of density n„/ 
driven by the ac electric field. Although these quasipar- 
ticles have their own kinetic inductivity,"'' their reactive 
contribution to the total kinetic impedivity is negligible 
at the frequencies of interest here. The quasiparticles' 
only significant contribution to the ac normal-fluid cur- 
rent density is therefore 

jni = (^nfE, (75) 

where 

(^nf = = n„/e^T„//m. (76) 

We can express the normal-fluid response in terms of the 
normal- fluid's kinetic impedivity Z^n = Pnf, so that E = 

^knjnl — Pnfjnf- 

Because the total ac current density carried by the 
strip is ji = jsi + jni, and E = Z^jsi, the overall kinetic 
impedivity Zk = TZk + i^k = 'R-k + of the strip, 

including the normal-fluid response, is the impedances- 
in-parallel combination, 

Zk = {Z^l+Z^^)-\ (77) 

At low temperatures the normal-fluid density is much 
smaller than the superfluid density, and \Zkn\ ^ \Zks\i 
such that Zk is dominated by the superfluid kinetic 
impedance. However, as T — > T^q, \Zks\ diverges, t„/ — > 
r, pnf Pn = l/cn, and Zk approaches the normal- 
state resistivity. 

Our results reduce to the well-known two-fluid 
description^* in the limit as jso — ^ 0, when the impediv- 
ity can be expressed in terms of the complex conductivity 
cr = CTi — i(T2, the linear response function connecting j 
and E calculated by Mattis and Bardeen.^^ As jsvi 0, 

Ti'ks 0, Xks — >• UjCks = MO'^'^O = l/'72, Cfnf — >■ CTi, 

Zka i/o-2-: l/cn, and {cri-i<J2)~^ = cr~^- 

At temperatures above T^q, the normal-state 
impedance is Z = i?„ -|- ioj{Lk + Lm), where 



is the geometric inductance associated with stored mag- 
netic energy and, for a long strip of conduction-electron 
density ric, total length £, width W, and thickness 
d, Lk = {m/nce^){i/Wd) is the normal-state kinetic 
inductance.^ The impedance is usually dominated by 
the normal-state resistance i?„ = pni/Wd except at 
very high frequencies.^ 



VII. DISCUSSION 

In this paper wc have presented fundamental theoreti- 
cal calculations of the kinetic impedance of thin and nar- 
row impure superconducting fllms for all temperatures 
and for all currents up to the depairing current. Our re- 
sults should be applicable to ongoing experimental stud- 
ies of small-scale superconducting devices in which the 
kinetic inductance plays an important role. 

Our results in the GL regime for the bias-current de- 
pendence of the kinetic inductance are in agreement with 
those of Anlage et al}^ and Annunziata et al? for the 
slow-experiment case and with the result of Anlage et 
al}^ for the fast-experiment case. However, Annunziata 
et al.,"^ in examining the case of T = and noting cor- 
rectly that Aq(0) is inversely proportional to Ao(0) (in 
our notation) [see Eq. (B3)], assumed that Aq(0) is in- 
versely proportional to Ag(0). However, this assumption 
is incorrect, as can be seen from Eqs. (23), (26), and (28). 
As a consequence, their prediction for the current depen- 
dence of the kinetic inductance does not agree with our 
results for either slow or fast experiments. 

In thin and narrow strips with sharp corners, ciirrcnt 
crowding leads to suppression of the order parameter in 
the immediate vicinity of sharp inner corners, and this 
can cause the critical current in such devices to be con- 
siderably lower than the depairing value. By op- 
timally rounding the inner corners, one should be able to 
raise the critical current to values close to the depairing 
value.19'23.51 

As discussed in Sec. II, our calculation of the criti- 
cal depairing current density has been carried out 
within a mean-field approach disregarding fluctuations. 
However, the experimental critical current density jc 
could be somewhat smaller than as a result of ther- 
mal or quantum fluctuations, which can initiate phase 
slips in ID wires^^'^'' or vortex nucleation at the edges 
of strips"'^^'"'^^'^'^'^^ when an energy barrier is overcome or 
suppressed to zero. The fluctuation-limited jc therefore 
may prevent the observation of both the predicted di- 
vergence of the slow-experiment kinetic inductivity at jd 
and the approach to the maximum values of the fast- 
experiment kinetic inductivity (shown by the filled sym- 
bols in Fig. 8). In fact, previous experimental observa- 
tions of a relatively small kinetic inductance rising to a 
peak and then rapidly dropping to smaller values with in- 
creasing applied current^'^*'^^'^'^ are explainable in terms 
of the growth of high-resistance normal regions as the 
current density rises above jc- 
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As discussed in Sees. III-V, relaxation of the supcrfluid, 
characterized by the relaxation time t^, plays an impor- 
tant role in determining both the real and imaginary 
parts of the kinetic impedance. To examine the physics 
of relaxation dynamics is beyond the scope of this paper, 
and for further discussion we refer the reader to Refs. 24 
and 29 and references therein. Nevertheless, Figs. 9 and 
10 suggest means by which Ts could be obtained from ex- 
perimental determinations of the superfluid's ac resistiv- 
ity and inductivity. At the very least, such experiments 
should be able to reveal whether they are in the slow- 
or fast-experiment limit, and experiments carried out at 
different frequencies and temperatures might be able to 
show the transition between these two limits. 

In this paper we have used the concept of the q- 
dependent penetration depth, which increases as the cur- 
rent density increases. This occurs because an increase 
of the current density causes a decrease in the magni- 
tude of the superconducting order parameter. This con- 
cept is the basis of the effective field-dependent pene- 
tration depth in type-II superconductors X{B,T), which 
diverges as B ^ Bc2{T). This quantity has been in- 
troduced to understand such phenomena as small-angle 
neutron-scattering form factors,^* magnetic coupling of 
vortex lattices in dc superconducting transformers, 
magnetization curves in low-pinning superconductors,^""^ 
elastic properties of the vortex lattice,®^ and ^SR mea- 
surements in the mixed state. 

Here we have used the wavevector q = 2'kAs/4)o to 
discuss the nonlinear response of the supercurrent to the 
gauge- invariant vector potential Ag in the local limit. It 
is important not to confuse our q with the q used by 
Tinkham^^ in his discussion of the nonlocal electrody- 
namics of superconductors, where it arises from Fourier 
transforming the convolution integral relating the linear 
response of the current density j{r) to the vector poten- 
tial A{r'), as calculated by Mattis and Bardeen.^^ 
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Appendix A: Unq 

The function Unq{ri, e, C) is the solution of the quartic 
equation 

e^u^ - 2T]eu^ + {if + e^- e^C^)u^ - 2??ew - r?^ = 0, (Al) 



obtained from Eq. (9). The desired solution of the quartic 
equation is 



Unq = {f + g + ??/e)/2, 



(A2) 



where 



a = -2[e6(c2- 1)3- 36^(1 + 16C2 + CW 

+ (A3) 

(A4) 
(A5) 
(A6) 



b = \/-4c'5 -f a2, 
c = e^(C^-l)-,;^ 
d = (a + 6)^3, 

22/3c2 + 



2^/ne^d 



(A7) 
(A8) 



^ 2[2e^(C^ 1) + ,^]_^^2(C^ + 1),_ ^^^^ 



3e2 



A similar solution of the quartic equation that arises 

in the closely related problem of the density of states in 
the Abrikosov-Gor'kov theory3^ was obtained in Ref. 65. 



Appendix B: Zero-current limit 

In the limit of zero current, g ^ (fg — >■ 0), m„o = 
and we obtain from Eq. (11) after defining t = T/Tcq and 
5q{T) = Ao(T)/27rfcsTeo, 



In 



1 °^ 1 

- = E(— 



n=0 



+ 1/2 ./{n+lJW+WW 



,(B1) 



which yields the temperature dependence of Ao(T) for 

all temperatures between and Tco, the transition tem- 
perature for (7 = 0. An analysis of this equation as i — )■ 
reveals that 47r5o(0) = 2i:^o{<S) /ksTcQ = 2Tie-'^ = 3.528, 
which is consistent with the fact that the above form 
of the Usadel theory coincides with the weak-coupling 
BCS theory for an s-wavc isotropic gap on a spherical 
Fermi surface. Values of Ao(T)/Ao(0) vs t = T/Tm, 
which reproduce well-known results, 33'66 can be obtained 
by numerically carrying out the sum in Eq. (Bl), and 
it can be shown analytically that [Ao(T)/Ao(0)]2 — >■ 
[8e2T/7C(3)](l -t) = 3.016(1 - i) as t -s- 1. 

In the zcTO-q limit, the sum in Eq. (16) can be evalu- 
ated analytically as shown in Eq. (14), such that'^^'^'^ 



nsojT) ^ A|(0^ ^ Ao(r) 
n,o(0) A2(T) Ao(0) 



tanh 



2kBT 



(B2) 



Prom Eq. (2), = Tiv p / it , and the normal-state 
conductivity C7„ = 2e^N{Q)D, we obtain 



nso{T)e^ 7ra„Ao(T) 



MoA§(T) 



m 



n 



tanh 



2kBT 



(B3) 
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t = TITco 

FIG. 11: Upper curve: reduced superfluid density 
n,o(T)/n,o(0) = [A(,(0)/A(,(T)]" vs t = T/T^o, obtained 
from Eq. (14). The dashed line, 2.660(1 - t), shows the 
slope as t — ^ 1. Lower curve: reduced upper critical field 
Bc2(T)/B,2(0) = [e(0)/e(T)]2 vs t = T/Tco, obtained from 
Eqs. (C9) and (CIO). The dotted line, 1.444(1 -f), shows the 
slope as t ^ 1. 




q/qm(0) 



FIG. 12: tcq = Tcq/Tco (solid) vs q/qm{0) obtained from Eqs. 
(23)-(25). Also shown are expansions of tcq about q = 
[dashed, Eq. (C3)] and q = qm{0) [dotted, Eq. (04)]. This 
figure is the same as a plot of t = T/Tco along the ordinate vs 
<?m(T)/iJm(0) along the abscissa. Aq{T) > only for values 
of t and g/gm(0) under the curve. 



Figure 11 exhibits the temperature dependence ofusoiT). 
As t ^ 1, n,o(r)/n.o(0) ^ [47reV7C(3)](l - t) = 
2.660(1 - 1). 
At T = 0, we can write 



Mo A, 



but because the zero-temperature London penetration 
depth Ai(0) can be expressed as^^ 



(B5) 



we see that Ao(0) = Al(0)(^o/^)^/^- (Recall that f < 
in the dirty limit under consideration here.) 



Appendix C: Zero-gap limit 



To find the boundary in the t-q plane where Aq(T) is 



reduced to zero, note from Eq. (9) that ey 1 + u'^q 

n + l/2 + Q/2TTkBT in the limit e = A,(T)/27rfcBT -> 0, 
such that Eq. (11) then yields the g-dependcnt tran- 
sition temperature. Defining tcq = Tcq/Tco, P = 
Q/27rkBTco, Pm = Q,n/2TTkBTco = e'/'(i/2) = e-yA = 
0.140, = nkBTcoe-y2 = hDql{0)/2, q^{0) = 

{TTkBTcoe-yhDy/^ = (3Mo^)'/', and Vm{0) = 
{h/2m){-nkBTcoe-'^ /hDy/'^, we obtain 



In 



1 



E 

n=0 



1 



n+ 1/2 



1/2 + P/tc 



t. 



cq 



, (CI) 



(C2) 



where ip is the digamma function and P/Pm = 
[9/9™ (0)]^ = [ws/fm(0)]^. Figure 12 shows tcq as a func- 
tion of q/qm{^) = 'Vs/vjn{Q). Expansions of tcq about 
5 = and q = qm{0) yield, respectively, the approxima- 
tions 



tc 



1 - 0.693 



(C3) 



9m(0) 



= 0.972 Jl- 



9m(0) 



,(C4) 



which are shown as the dashed and dotted curves in Fig. 
12. 

A similar procedure can be used to determine qm{T), 
the value of q that drives Ag(T) to zero for a given value 
of t. Equation (11) then yields 



In- 
t 



^ 00 



?i=0 



n + 1/2 n+l/2-h( 



= ^ 



(C5) 
(C6) 



where a is given in Eq. (34). Since tcq and qm{T) are 
both determined by the equation obtained by setting 
Ag(T) = 0, it should not be surprising that a plot of 
t vs qm{T)/qm{0) is exactly the same as the plot of tcq vs 
<l/<lm{Q), shown in Fig. 12. Expansions of qm{T)/qm{0) 
about t = and t = 1 yield, respectively, the approxima- 
tions 



5m (0) 



e27 

= 1 - —t"^ = 1 - 1.057t2 



(C7) 




e-'^VT^ = 1.202vT^, (C8) 



which correspond to the dotted and dashed curves in Fig. 
12. Equations (C7) and (C8) are most easily obtained by 
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making the replacements q — )• qm{T) and tcq t in Eqs. 
(C4) and (C3). 

The upper critical field is related to the temperature- 
dependent coherence length via Bc2{T) = ^o/27r^^(T), 
and in the dirty limit Bc2 can be obtained from the 
equation' 



68 



In 



where t = T/Tco and 



(?)( 



(C9) 



(CIO) 



where 5^2(0) = (f>o/2TT£,^{0) and ^(0) = (7r^o^/3)^/^ 
At t = 0, we have Ao(0)/C(0) = vW^^i(O)/^ = 
0.977Al(0)/£. 

Ast^O, B,2{T)/B,2{0) « 1 - (2eV3)t2. 

As t ^ 1, B,2{T)/B,2{0) ^ (8e'^/7r2)(l - = 
1.444(1 — f). Note that this result is consistent with what 
Helfand and Werthamer^^ found for their normalized 
field h*{0) = Ba2{0)/{dBc2/dt)t=i = 0.69 = 1/1.444 in 
the dirty limit. Since ^(T) 



^7r3e-T/24^/^/^/^^ = 
0.852VCo?/Vl - t, we have in this limit Ao(T)/^(T) = 
K = V42C(3)/7r4AL(0)/£ = 0.720Al(0)/1 

A plot of Bc2{T) / Bc2{0) is equivalent to a plot of 
K(0)/C(T)]2, as shown in Fig. 11. Moreover, compar- 
ing Eqs. (C9) and (C9) with Eqs. (C5)-(34), we see that 
qm{T) = 1/C(T) and ^^(O) = 1/^(0), and for all temper- 
atures we have 



Bc2{0) e{T) qliO) ■ 



(Cll) 



Appendix D: Work done and free-energy changes 
resulting from current changes 

It is of interest to examine the changes in energy as 
the current increases from zero to some final value. The 
work done per unit volume is 



j,g,Edt' = -(^) f'jsg'dq', 



(Dl) 



which is equal to the change in the free-energy density, 
as we show below. 

The free-energy density 17g(T) = Fs{T) - Fn{T) of 
a current-carrying superconductor relative to the energy 
density of the normal state can be obtained by taking 
advantage of the theoretical similarities to the problem 
of superconducting alloys containing paramagnetic im- 
purities. The expression for f],(T) = Fs{T) - Fn{T) for 
the latter case obtained by Skalski et al.^° in their Eq. 
(5.6) can be rewritten compactly for the current-carrying 
superconductor as 



Slq = N{0) [ Re 2p(w)-|- 



tanh^dw, (D2) 



using the replacements and changes in notation F — > Q, 
A(T,r) Aq{T), No N{0), and ^ oo (weak- 
coupling limit). Here /3 = l/ksT, 



= u 1 



c 



and 



1 - i 



doj' 



c 



l-u) 



c 



2(w2_l) 



(D3) 



(D4) 



(D5) 



is a quantity arising from partial integration of the term 

proportional to ln(l + e~^") in the entropy contribution 
to riqiT). Some useful relations arc Re[p(0)] = and 



-(«2_ 1)1/2 VJ- 2- 

The quantity within the brackets in the integrand of 
Eq. (D2) can be expressed in terms of u and u as 



1 



1 



(D7) 



and the integral in Eq. (D2) can be evaluated using con- 
tour integration around the boundaries of the first quad- 
rant of the complex uj plane, taking into account the poles 
along the imaginary axis at the Matsubara frequencies 
icOn = i2-KkBT{n -|- 1/2) = i^nkBTr]. The result is 

n,(r) = -A^(0)(27rfcBT)^f;J2( - l) 



+ - 



(D8) 



Diflferention of Eqs. (11) and (D8) with respect to q 
yields the general result that 



dq 



= N{0){2'KkBT) 



dQ 
dq 



CO ^ 



(D9) 
(DIG) 



This result agrees with Eq. (Dl). 

In the zero-temperature limit, the sum in Eq. (D8) can 
be converted to an integral over 77 using Eq. (9) , with the 
result^^ 



2 

N{0)Aq{0r 



1 







2 y'-^^ 



2 3 
3 



, Co < 1, (Dll) 



(2Cg + l)v^ 
3Co 



+Cotan-\/Co2-l 



Co>l, (D12) 
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FIG. 

Eqs. 



13: Energy densities at T = 0: 0,(0) 
(Dll) and (D12)], Tcq{0), [dotted, Eq. 



[solid curve, 
(D13)], and 

TkgiO) [dashed, Eq. (D14)], normalized to N{0)Al{0)/2. The 
black points identify the values of these energy densities at 
9d/9m(0) = 0.689, where f2g(0) has its maximum slope and 
the depairing current density is achieved. 



where Aq(0) is given by Eqs. (23) and (24), and Co by 
Eq. (25). 

It is possible to write the free-energy density as the 
sum of two terms, flq{T) = Fcq{T) + Tkg{T), where we 
identify TcqiT) as the condensation-energy density and 
J^kq{T) as the kinetic-energy density: 

oo 

jr,g(T) = -iV(0)(27rfcBT) V hhcoj /"^ : - l) 

\/l + <a ^ 



+ 



1 



(D13) 



OO ^ 

J^kq{T) = N{0){2nkBT)Q ^ — (D14) 
Note, however, that both terms depend upon q and inter- 



act. As q increases from zero, Tkq{T) initially increases 
from zero and the magnitude of Tcq{T) decreases, but 
their sum ^q{T) decreases to zero as g ^ qm{T). Figure 
13 shows the q dependence of the energy densities Q.q, 
Fcq, and Fkq at T = 0. 

It is important to note that Tkq{T) is not equal to 
Ck{q-iT)j1 12 except in the limit bs q^Q [see Eq. (1)]. 



Appendix E: Bulk thermodynamic critical field 

The bulk thermodynamic critical field Hc{T) is defined 
via the superconducting condensation energy at g = 0: 



1 



(El) 



= iV(0)(27rfcBT)2^(2 



n=0 



-2r?- 



where t? = n -|- 1/2 and cq = Ao(r)/27rfcBr. Equa- 
tion (E2), which follows from Eq. (D8) when (j — > and 
Unq — > Uno = rj/co [sce Eq. (9)], reproduces the BCS'^^ 
temperature dependence of Hc{T). When T ^ 0, the 
sum can be replaced by an integral over rj, and the result 
is 



Ii^oHHo) = ^Ar(0)A^(0), 



as expected from Eq. (Dll). As t = T/Tco 1, 



(E3) 



Hc{T) 



HM L7C(3) 



1/2 



{l-t) = 1.737(1 - t). (E4) 
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